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ABSTRACT 

We study the ionization and kinematics of the ionized gas in the nuclear region of the barred Seyfert 2 galaxy NGC 5643 using MUSE 
integral field observations in the framework of the MAGNUM (Measuring Active Galactic Nuclei Under MUSE Microscope) survey. 
The data were used to identify regions with different ionization conditions and to map the gas density and the dust extinction. We 
find evidence for a double sided ionization cone, possibly collimated by a dusty structure surrounding the nucleus. At the center of 
the ionization cone, outflowing ionized gas is revealed as a blueshifted, asymmetric wing of the [OIII] emission line, up to projected 

velocity v 10 -450 bn s -1 . The outflow is also seen as a diffuse, low luminosity radio and X-ray jet, with similar extension. The 

outflowing material points in the direction of two clumps characterized by prominent line emission with spectra typical of HII regions, 
located at the edge of the dust lane of the bar. We propose that the star formation in the clumps is due to ‘positive feedback’ induced 
by gas compression by the nuclear outflow, providing the first candidate for outflow induced star formation in a Seyfert-like radio 
quiet AGN. This suggests that positive feedback may be a relevant mechanism in shaping the black hole-host galaxy coevolution. 

Key words. Galaxies: active - Galaxies: individual (NGC5643) - ISM: jets and outflows - Techniques: imaging spectroscopy 


1. Introduction 

Active Galactic Nuclei (AGN) have a profound influence on their 
host galaxies, as they are capable of ionising large fractions of 
the interstellar medium (ISM) and of accelerating fast outflows. 
These are powerful enough both to sweep away most of the gas 
in the ISM and to heat the gas in the surrounding halo, inhibiting 
gas accretion (“negative feedback”, see Fabian et al. |2012l and 
references therein). Fast outflows could also induce star forma¬ 
tion (“positive feedback”) in the pressure compressed molecular 
clouds (e.g. Silk et al. 120131 or in the outflowing gas (Ishibashi 
& Fabian l20T2l20T4l Zubovas & King l20l4l ). 


* This work is based on observations made at the European Southern 
Observatory, Paranal, Chile (ESO program 60.A-9339) 


Although such ‘positive feedback’ has been invoked in re¬ 
cent years by several theoretical works to explain observed cor¬ 
relations between, e.g., AGN luminosities and nuclear star for¬ 
mation rates (Imanishi et al. [20111 Zinn et al. 120131 Zubovas et 
al. 120131) . black hole accretion and star formation rates in AGN 
(Silverman et al. 120091 Mullaney et al.|2012 Silk l2013l) or down¬ 
sizing of both black holes and spheroids (Silk & Norman (20091 . 
still very few observations of such phenomenon at work have 
been argued for, mostly related with spatial alignment of star 
forming regions or companionion galaxies with very powerful 
radio jets (e.g. Kramer et al. 120041 Croft et al. 120061 Elbaz et al. 
2009, Feain et al. 120071 Crockett et al. 120121 Salome, Salome 
& Combes 12015) . The only known example of “positive” feed¬ 
back at high-z in a radio quiet AGN has been recently presented 
by Cresci et al. (120151 ). who detected with near-IR Integral Field 
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Fig. 1 . The galaxy NGC 5643. Panel a: Three color (Blue: B band; Green: V band; Red: I band) image of the galaxy obtained with the 2.5m du 
Pont Telescope for the Carnegie-Irvine Galaxy Survey (CGS, Ho et al. l201 it . The MUSE field of view is shown as a white square. North is up and 
East is left. Panel b: Three colors MUSE map, showing Hff emission in Blue, [OIII]/15007 emission in Red and the integrated continuum emission 
from 5200 to 6000 A in Green. The tip of the dust lane in the bar, visible on larger scale in Panel a, is evident next to the blue clump East of the 
nucleus as a deep in the continuum. Panel c: Three colors MUSE maps, showing Hu in Blue and [OIII]/15007 in Red as before, and [NII]T6584 
emission in Green. The two-sided AGN ionization cone traced by the filamentary [OIII] emission is prominent, as well as the strongly Ho: emitting 
clumps between the two eastern lobes of the cone. Nord is up and East is left. 


observations star forming clumps in the host galaxy of a z ~ 1.6 
quasar, possibly triggered by the outflow pressure at its edges. 
However, we are still lacking observations of positive feedback 
at lower AGN luminosity and outflow energy, in order to under¬ 
stand its importance in normal AGNs. 

In this paper we present observations of the nuclear region 
of NGC 5643, a barred, radio-quiet Seyfert 2 galaxy at a dis¬ 
tance of ~ 17.3 Mpc seen almost face on (i ~ -27 ± 5 deg, de 
Vaucouleurs et al. 11976 1 . with a well known ionization cone ex¬ 
tending eastward of the nucleus, parallel to the bar (Schmitt et 
al. 1 19941 Simpson et al. 1 19971 Fischer et al. l20131 l. A three color 
image (B,V, I bands) from the Camegie-Irvine Galaxy Survey 
(CGS, Ho et al. 2011 ) is shown in Fig. [Tj panel a. As usually 
observed in barred spiral galaxies, sharp, straight dust lanes ex¬ 
tend from the sides of the central nucleus out to the end of the 
bar, roughly parallel to its major axis and being more clearly vis¬ 
ible to the East of the nucleus. The galaxy is also known to host 
a diffuse, low luminosity radio jet on both sides of the nucleus, 
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which is nearly 30” long (2.5 kpc , Leipski et al. 120061 Morris et 
al. fl985) . 

This galaxy was observed with the Multi Unit Spectroscopic 
Explorer MUSE, the optical large field integral-field spectrom¬ 
eter at the VLT (Bacon et al. 120101 as part of the science ver¬ 
ification run, in the framework of the MAGNUM (Measuring 
Active Galactic Nuclei Under MUSE Microscope) survey, a pro¬ 
gram aimed at the observation of nearby Active Galactic Nuclei 
(AGNs), to study the physical conditions of the Narrow Line 
Regions (NLRs), the interplay between nuclear activity and star 
formation, and the effects and acceleration mechanisms of out¬ 
flows. 

In the following we adopt the systemic velocity measured 
by Koribalski et al. (120041 . 1199 km s 1 . We use Hq = 
69.6 km s ! Mpc 1 and O m = 0.286 (Hinshaw et al. 120131) . 
At the distance of the galaxy of 17.3 Mpc, the angular scale is 
1” = 83 pc. 
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Fig. 2. Continuum subtracted, rest frame spectra extracted from a region of 1.6” x 1.6” around the nucleus (upper panel) and around the star 
forming clump A {lower panel). The different excitation between the two regions is evident from the different line ratios, where [OlII] and [Nil] 
are brighter in the nucleus, while stronger Her and H/3 emission is detected in the star forming clumps. 


2. Observations, data reduction and analysis 

NGC 5643 was observed with MUSE during the science ver¬ 
ification run on June 24, 2014, under program 60.A-9339 (PI 
Marconi/Hawthorn/Salvato). The average seeing during the ob¬ 
servations, derived directly from foreground stars in the final dat- 
acube, was FWHM = 0.88” + 0.02”. The nuclear region of the 
galaxy was observed with four dithered pointings of 500 s each, 
with the sky sampled with 4 shorter 100 s exposures in between. 
The data reduction was performed using the recipes from the 
early-release MUSE pipeline (version 0.18.1), as well as a col¬ 
lection of custom IDL codes developed to improve the sky sub¬ 
traction, response curve and flux calibration of the data. Further 
details on the data reduction will be provided in a forthcoming 
paper (Carniani et al., in preparation). The final datacube con¬ 
sists of 315 x 315 spaxels, for a total of almost 100000 spectra 
with a spatial sampling of 0.2” x 0.2” and a spectral resolution 
going from 1750 at 465 nm to 3750 at 930 nm. 

Moreover, NGC 5643 was observed in the X-rays by the 
Chandra Observatory on December 26, 2004, for 7.5 ks with 
the Advanced CCD Imaging Spectrometer (ACIS: Garmire et al. 
120031 . Data were reduced with the Chandra Interactive Analysis 
of Observations 4.5 (CIAO; Fruscione et al. 120061 and the Chan¬ 
dra Calibration Data Base 4.6.5, adopting standard procedures. 
In order to increase the spatial resolution, the imaging analysis 
was performed by applying a factor of 2 sub-pixel event reposi¬ 
tioning (e.g. Wang et al. 12011S >. We therefore used a pixel size of 
0.246”, instead of the native 0.495”. 


2.1. Modeling the stellar continuum and emission lines 

In presence of young/intermediate-age stellar populations, which 
present intense Balmer absorption (e.g. EW(H/3) can reach up to 
8A), a careful subtraction of the underlying stellar continuum is 
important in order to recover the correct emission line fluxes and 
profiles. Although the optimal approach is to make a simulta¬ 
neous fit of the stellar continuum and of the emission lines (see 
e.g. Sarzi et al. 12006 1). with 10 5 spaxels and a spectral sampling 
of 1.25A this becomes computationally unfeasible for a MUSE 
datacube on a standard desktop computer. On the other hand, 
considering the intensity of the emission lines we aim at study¬ 
ing here (see next Section), sufficient accuracy can be obtained 
by fitting the stellar continuum while avoiding the regions in¬ 
terested by emission lines and then by extrapolating the fit over 
those regions. We use the pPXF code (Cappellari & Emsellem 
12004) to perform this task over the wavelength range of interest, 
i.e. shortward of 7000A. A linear non-negative combination of 
stellar spectral templates is fitted to the data by adjusting the sys¬ 
temic velocity and the velocity dispersion together with the co¬ 
efficients of the linear combination. We consider 40 stellar tem¬ 
plates from the simple stellar populations provided by Bruzual 
& Chariot (120031 0. distributed in a grid of 10 ages (5.2, 25, 100, 
290, 640, 900 Myr, 1.4, 2.5, 5.0 and 11 Gyr) and 4 metallic- 
ities (0.2, 0.4, 1, and 2.5 times solar). Each spectral template 
is converted from wavelength to velocity space and so is each 
observed spectrum. pPXF iteratively adjusts the systemic veloc- 
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ity and the gaussian velocity dispersion of each template (this is 
done by plain convolution in the velocity space with a gaussian 
kernel; velocity and dispersion are assumed to be common to all 
component templates) and computes the coefficients of the lin¬ 
ear combination of templates that best reproduce the observed 
spectrum. The procedure is iterated until the minimum \ 2 and 
the best fitting systemic velocity, velocity dispersion and coeffi¬ 
cients of the template linear combination are obtained (see Cap- 
pellari & Emsellem 120041 for more details). Note that the tem¬ 
plate distribution in age and metallicity allows to cover virtually 
any (combination of) stellar populations inside a galaxy, includ¬ 
ing the very young ages that characterise the star-froming blobs. 
We mask regions corresponding to ±750km s 1 around the main 
astrophysical emission lines and around the Nal absorption at 
rest and in the galaxy’s restframe to avoid contamination from 
the interstellar absorption. We note that the fit is expected to be 
definitely degenerate in terms of the possible linear combina¬ 
tions of SSPs, but this is irrelevant to our goal: a realistic model 
of the star formation history and metallicity distribution is not re¬ 
quired to ensure just an accurate stellar continuum subtraction. 
We verify that we obtain an accurate subtrction by comparing the 
residuals from the fit with the estimated noise level and finding 
them consistent with each other. 

This fitting procedure produces reliable results only if a suf¬ 
ficient signal-to-noise ratio (SNR) in the continuum is provided, 
namely > 10 per A. In the original datacube at full resolution 
this is achieved only in the brightest parts; for the fainter regions 
we apply an adaptive smoothing procedure which extends the 
adaptsmooth code (Zibetti [2009al Zibetti et al. 12009b) to work 
with spectral datacubes (azmooth3, Zibetti et al. in preparation). 
At each spaxel the SNR is estimated in a narrow wavelength win¬ 
dow on a featureless continuum, between 6000 and 6050 A. If 
the SNR is lower than the threshold of 20 per spectral pixel, con¬ 
centric annuli of spaxels are coadded iteratively until the SNR of 
20 is reached. No more than 10 annuli are coadded at each po¬ 
sition. The smoothed cube warrants sufficient SNR at all spaxels 
while preserving as much as possible the effective spatial resolu¬ 
tion. The typical smoothing radius (in pixels, 1 pix=0.2”) varies 
from 1 (i.e. no smoothing) within 5” from the centre to 2, 3, and 5 
at distances of ~ 13”, 18”, and 25”, respectively. Larger smooth¬ 
ing radii are required only in sparse regions, typically outside 
25”. Indeed, most of the analysis presented in this work is unaf¬ 
fected or marginally affected by the smoothing: the EELR are in 
fact mostly contained in the no-smoothing zone or in the 2-pixel 
radius smoothing zone. In any case, the smoothing affects scales 
typically < 2” in diameter, hence smaller than those we investi¬ 
gate in the following. The output is then processed with pPXF as 
explained above to produce a pure stellar continuum cube, stellar 
velocity and stellar velocity dispersion maps. 

2.2. Emission line fitting , density and extinction maps 

The obtained stellar continuum is subtracted from the original 
data, and the resulting datacube is used for the following analy¬ 
sis. We fit the main emission lines accessible in the MUSE wave¬ 
length range, i.e. H/3 , [OIII]T4959, [OIIIJT5007, [OIJ/16300, 
[NII]/16548, Her, [Nil]46584, [SIIJ46716, [SII]46730. To better 
reproduce the observed spectral profiles in the data, we use both 
a two Gaussian components fit and a single Gaussian fit to each 
line. The flux obtained with the double Gaussian is used in the 
following only for the spaxels where the second component has a 
SNR > 3 in the fit, i.e. in the nuclear region where the line pro¬ 
file is more complex (see Fig. [2}. This choice allows us to use 
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the more degenerate two components fit only where it is really 
needed. 

From these fluxes we constructed emission line maps in the 
Ha [OIII] and [Nil] transitions (see Fig. [TJ,c). These maps are 
already clearly showing that the line emitting gas presents an 
elongated morphology, aligned with the bar in the E-W direc¬ 
tion. We will refer in the following to this region as Extended 
Emission Line Region (EELR). As reported by Schmitt et al. 
(11994b . the EELR is elongated both to the E and to the W of 
the nucleus, although it is much fainter to the W, especially at 
shorter wavelengths (e.g. Fig. [T] panels, b and c). As we will dis¬ 
cuss later, this is due to the higher extinction to the W, indicating 
that the EELR is tilted with respect to the galaxy plane, and it 
is lying behind the disk of the galaxy to the W, while it is seen 
mostly above the disk to the E. 

An [OIII] emitting structure is very prominent, shown in red 
in Fig. [] originating at the AGN location in the galaxy center 
and extending in the E-W projected directions, part of an ion¬ 
ization bi-cone (see Fig. |T] panels b and c; see also Fischer et 
al. 120131 ). A V-shaped ionization cone in the core of NGC 5643 
was already reported by Simpson et al. ( 119971 ) using HST nar¬ 
rowband imaging, although the western part of the bi-cone was 
not detected in their data, probably due to lower sensitivity. 

Strikingly, a blue. Her bright region is located exactly at the 
tip of the Eastern part of the EELR at a distance of ~ 15” = 
1.2 kpc from the nucleus, at the location where the line emitting 
region reaches the dust lane of the bar (clump A, see panel b, 
where the dust in the bar is visible as a dip in the green contin¬ 
uum), and another Ha bright clump is located ~ 8” to the south¬ 
west of clump A (clump B, see also Fig.[3]and Fig. [7]). A compar¬ 
ison between the line ratios measured on clump A and the line 
ratio in the nuclear region shows that the excitation mechanisms 
are different, with the nucleus showing an AGN-like ionization 
as expected, while clump A displays the typical spectrum of a 
star forming region (see Fig. 2]). We will discuss further these 
star forming regions in the following. 

We estimate the electron density from the ratio between the 
[SII]/16716 and [SII]46730 lines, which is a widely used in¬ 
dicator of this quantity (e.g. Osterbrock & Ferland [2006 1. We 
compute the line ratio in each spaxel, and convert it to an elec¬ 
tron density using the IRAF task temden, assuming a temper¬ 
ature of 10 4 K. The resulting electron density map is shown 
in Fig. [3]a, for the regions where the [SII] lines are detected 
with SNR > 5. The density is higher in the nucleus, and it de¬ 
creases outwards by almost three order of magnitudes to a level 
of ~ 10 cm 3 . The density is again higher (~ 100 an 3 ) in the 
blue, Ha bright clump. We note that the density derived with this 
method on the nucleus is higher than the critical density for [SII] 
(~ 3 x 10 3 cm 3 ) above which the lines become collisionally de- 
excited (see Osterbrock & Ferland[2006). Therefore, the density 
derived in the nucleus with this method is very uncertain. 

We also use the Balmer decrement Ha/H/3 to derive the dust 
extinction map, assuming a Calzetti et al. (120001 ) attenuation law 
and a fixed temperature of 10 4 K. Although there could be in 
principle a large temperature range across the field of view, the 
expected variation in Ay is +0.2 for a reasonable temperature 
range between 5000 K to 20000 K, typical of the conditions in 
HII regions and in the Narrow and Broad Line Regions of AGNs 
(Osterbrock & Ferland 120061 ). This do not change our qualitative 
picture as the difference between the dusty structure W of the 
nucleus and the EELR to the E is one order of magnitude greater, 
the expected variation for a temperature of 5000 K or 20000 K 
is ±0.2 in Ay. which do not change our qualitative picture as the 
difference between the dusty structure W of the nucleus and the 
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Fig. 3. Electron density (panel a) and extinction map (panel b), as de¬ 
rived from the [SII]/16716/[SII]/16730 and Hcr/HyS line ratios (see text) 
with Ho- contours overplotted. The star forming clump has higher den¬ 
sity and extinction than its surroundings. A dusty structure with very 
high extinction is detected just West of the nucleus and Her peak, corre¬ 
sponding to the tip of the bar driven dust lane. Both maps are displaying 
the spaxels with SNR > 5. The ratio between the continuum flux in the 
8000-8500 A spectral region and the 5100-5500 A one, with Her con¬ 
tours overplotted, is shown in panel c. While the star forming clumps 
and the nucleus are bluer due to the intrinsic shape of the continuum, 
the higher extinction regions are showing redder colors, tracing the dust 
distribution in the regions even where the SNR of the emission lines is 
lower. The dust lanes to the W and E of the nucleus, an extinction ridge 
to the S connecting the dust lane with the dust structure surrounding the 
nucleus are clearly visible. Panel d shows the rest frame Her Equivalent 
Width map. The star forming clump A. marked with a circle as well 
as clump B, shows the highest £VT(Hcr) = 350 A in the field of view, 
suggesting an age of ~ 10 7 yrs (see text for details). 


EELR to the E is one order of magnitude greater. The resulting 
map is shown in Fig. [3jb. We find that the reddening is generally 
higher to the W of the nucleus. This confirms that the bright line 
emitting gas is probably equally extended to the W and to the E 
of the nucleus, but while the E part is seen in front of the bar, 
the W side is hidden behind the bar material and suffers higher 
extinction. 

Moreover, the dust extinction is highest in a curved structure 
just W of the nucleus, which seems to be connected to the dust 
lane in the bar through an higher extinction ridge, visible at the 
bottom part of the map and around clump B (see also Fig. [3 c, 
and Fig. [T| panel b, where the dust lane is visible as a dip in the 
continuum shown in green). This dusty structure appears to be 
~ 7” x 1.5” in size, corresponding to 580 x 125 pc, and shows 
the highest reddening in the central region of the galaxy, with 
A v > 3 and a maximum of Ay = 4.8. The stmeture is also ev¬ 
ident in Fig. [T] as a dip in the shorter wavelength emission line 
map, [OIII]/15007 shown in red, as well in Fig. [3}:, where the 
ratio between the continuum flux in the 8000 - 8500 A spectral 
region and the 5100 - 5500 A one is shown. While the star form¬ 
ing clumps and the nucleus are bluer due to the intrinsic shape 
of the continuum, the regions with higher extinction have redder 
colors, tracing the dust distribution in the regions even where the 
SNR of the emission lines is lower. In this continuum ratio map 
the dust lanes to the W and to the E of the nucleus are clearly 
visible, as well as the higher extinction ridge to the S connecting 
the dust lane with the dust structure surrounding the nucleus. It 
has been suggested that this dusty region represents the tip of 
the bar in the nuclear region, where the flowing material forms 
a disky structure elongated perpendicularly to the bar itself (e.g. 
Morris et al. 119851) . This is actually predicted by dynamical mod¬ 
els, where the gas is expected to settle into a disk with rotation 
axis along the bar (see, e.g., Durisen et al. 119831 . It is possi¬ 
bly also responsible for the collimation of the ionizing photons 
along the bar (Menezes et al. 1 2015 ). generating the ionization 
bi-cone seen in [OIII] line emission. The alignment between the 
inner bar and the EELR is therefore actually expected, with the 
ionizing continuum that escapes in this case in a cone slightly 
tilted with respect to the bar so that the gas is mostly illuminated 
in front of the bar to the E and behind the bar to the W. This 
high extinction region is also dominated by shock ionization, as 
discussed in Sect.0 

3. Resolved BPT diagrams 

We use the measured emission line fluxes in each spaxel to com¬ 
pute excitation maps of the galaxy region within the MUSE field 
of view. We consider both the classical Baldwin et al. (1198ll) 
BPT diagram using the [OIH]T5007/H/? versus [Nil]T6584/Ha 
line ratios (N-BPT in the following), as well as the alternative 
versions using [SII]T6716,30 (S-BPT) or [OI]T6300 instead of 
[NII] /16584 (O-BPT, see e.g. Kewley et al. 120061 Lamareille 
12010) . In these diagrams, galaxies dominated by AGN (Seyfert- 
type). Low Ionization Nuclear Emission line Regions (LINERs), 
shocks and star formation ionization populate different regions. 
For the N-BPT diagram, we also plot the location of points for 
two prototypical galaxies, the shock excited region of NGC 1482 
and for the AGN excited emission of NGC 1365, from Sharp & 
Bland-Hawthorn (120101) . The bisector between these two fiducial 
traces is used to discriminate between AGN excited and shock 
excited gas in the N-BPT diagram. In the right panels of Fig. [4] 
we show the maps where different regions have been colour 
coded according to their line excitation as derived from the cor¬ 
responding BPT diagram on the left. HII regions are marked in 
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Fig. 4. Resolved BPT diagrams for the nuclear region of NGC 5643. The N-BPT ([Nil]46584/Ha vs [OIII],l5007/H/3, upper panels), S-BPT 
([SIIJ/16716,30/Ha vs [OIII]/15007/H/3, central panels) and O-BPT [OI]/16300/Ha vs [OIII]/15007/H/J, lower panels) diagrams for each spaxel 
with S NR > 3 in each line are shown on the left, along with contours showing the density of the points on the diagram. The location of the Ha 
bright clumps A and B, as derived from an integrated spectrum extracted on a region of 1.6 x 1.6”, is shown by an yellow circle, while a cross is 
marking the location of the nucleus. The dashed, cyan lines in the N-BPT mark the location of points for the shock excited region of NGC 1482 
and for the AGN excited emission of NGC 1365 from Sharp & Bland-Hawthorn (2 01 Oi l. The bisector between the two, marked with a dashed red 
line, divide AGN dominated by shock dominated regions. A map marking each spaxel with the color corresponding to the dominant excitation 
at its location is shown on the right, again for spaxels with SNR > 3. The star forming regions are marked in blue, the Seyfert-type ionization 
is shown in green, LINER/shock dominated regions are shown in red and intermediate regions are shown in magenta in the N-BPT. Star forming 
regions are prevalent in the outer part of the MUSE field of view, AGN ionization is dominating in the core and in the ionization cone, while shock 
excitated gas is prevalent around the EELR and the cone. The nucleus is marked with a cross, and the Ha bright clups A and B are marked with a 
yellow circle. Clump A is dominated by star formation in all the three diagrams, while Clump B is classified as intermediate by the N-BPT. 


blue, AGN dominated regions are in green, LINER-like/shock 
excited regions are marked in red, while composite regions in 
the N-BPT are shown in magenta (see Kewley et al. 120061 Kauff- 
mann et al. 120031 ). In all BPT maps the nuclear region, the EELR 
and the ionization cone are dominated by AGN ionization. Shock 
excited gas is prevalent in the surrounding region and at the lo¬ 
cation of the dusty structure surrounding the nucleus, while star 
formation is dominant at the edges of the MUSE field of view, 
in a circumnuclear star forming ring. The only notable exception 
to this schematic picture are the blue, Ha bright regions already 
prominent in Fig. [T| Although they are located at the tip of the 
central AGN dominated region and along the ionization cone, 
they are classified as star formation dominated by all the three 
BPT diagrams (clump A), or by two out of three (clump B, see 
Fig.[4]and Fig. [7j. Interestingly, comparing the location of clump 


A and B in Fig. |T| they appear to be also located on the ridge of 
the dust lane in the bar. 

Clump A is also the region with the highest Ha EW in teh 
field of view, EW(Ha) = 350 A rest frame (see Fig. [3]d).The 
other star forming regions, mostly located in a star forming 
ring near the edge of the MUSE field of view do not exceed 
EW(Ha) = 250 A. The high EW in Clump A suggests a young 
age of its stellar population. Using the Starburst99 models (Lei- 
therer et al. 11999k assuming a solar metallicity and a Salpeter 
IMF, we obtain an estimate of the age of the cluster of tduster = 
5.5 x 10 6 yr for an instantaneous burst, and t c i uster = 8.7 x 10 1 yr 
for a continuous burst of star formation. These values are to be 
considered as upper limits, as the measured continuum under 
the line is integrated over a line of sight crossing of the whole 
galaxy, not only of the star forming region, and therefore the 
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Wavelength (A) 


Fig. 5. Two components Gaussian fitting of the 1.6”x 1.6” rest frame nu¬ 
clear spectrum of NGC 5643 for the Hy3,[OIII]T4959 and [OIII]/15007 
emission lines. The [OIII]/15007 lines in the nucleus shows an asymmet¬ 
ric blue wing, reproduced with blueshifted Gaussian components. The 
systemic components have a FWHM sys = 286±6 km , while the out¬ 
flowing components have a centroid shifted by Vo, ou < = -200 km s _I and 
a FWHM oul = 590 km s -1 . The projected velocity at the 10th percentile 
of the overall emission-line profile is Vio = -430 km s -1 , confirming the 
presence of outflowing, high velocity gas from the nuclear region. 


measured EW(H<r) is to be considered as a lower limit to the 
true value. From the measured Her flux extracted in a region 
of 1.6” x 1.6” around the clump, assuming the star formation 
rate calibration of Kennicutt & Evans (120121 1. and correcting for 
dust extinction using the Balmer decrement H<r/Hj6 assuming a 
Calzetti ( |20001 > extinction law, we derive a star formation rate 
in the clump of SFR( c i ump) — 0.03 M Q /yr (corresponding to 
2-sfr = 1.7 M 0 yr~ x kpc~ 2 ). Given that clump B is classified as 
‘intermediate’ by the N-BPT, we do not measure a SFR from Her, 
as it may be contaminated by AGN ionization. The spectrum of 
the clumps is characterized by an almost featureless continuum 
as expected from the cluster age estimated from the EW(Hcr). 
One might possibly analyze the Balmer lines in absorption in or¬ 
der to derive a more reliable stellar age, but with such a strong 
infilling emission (EW(H<r)>300A) one would require a much 
higher SNR in the continuum than actually available in order to 
reliably measure the Balmer absorptions. Actually, the emission 
is so strong that completely obliterates any absorption in these 
clumps. We will further discuss in the next Section why these 
particular locations are so special, and what has probably trig¬ 
gered the star formation in the clumps. 


4. Outflowing gas and positive feedback 

The presence of high velocity gas in the EELR region is revealed 
by the ionized gas dynamics traced by our MUSE observations. 
The nuclear spectrum shown in Fig. [2] presents an asymmetric 
profile of the [OIII] lines. The [OIII] lines are an ideal tracer of 
extended outflowing ionized gas, as they cannot be produced in 
the high-density, sub-parsec scales typical of AGN Broad-Line 
Regions (BLR), and have been widely used in the literature to 
study the winds properties (e.g. Villar-Martm et al. l201ll Zhang 
et al. 120111 Brusa et al. 1 2015 1. A similar, asymmetric profile is 
observed for other forbidden lines, as the [Nil] and [SII] dou¬ 
blets. 

We fit the H/i, [OIII]T4959 and [OHI]/15007 nuclear spec¬ 
trum of NGC 5643 (in a region of 1.6” x 1.6”) with a dou¬ 
ble Gaussian for each line, a systemic component and an out¬ 
flowing one (see Fig. [5]). All the lines are fitted simultane¬ 
ously, fixing the wavelengths ratios and [OIII] flux ratios ac¬ 
cording to theoretical values. The systemic components have 



[OIII] blue wing contours 

b 



Radio contours 



Fig. 6. The outflowing gas in the core of NGC 5643. Panel a: Ha map 
of the central region, with superimposed the contours of the flux in the 
blue wing (v < -400 km s _1 ) of the [OIII]/t4959 emission line, trac¬ 
ing the high velocity outflowing ionized gas. Panel b: Ha map with 
superimposed the 8.4 GHz radio contours from the VLA observations 
presented in Leipski et al. 120061 The diffuse radio jet is coincident with 
the location of the blue [OIII],14959 wing. Panel c: Ha map with su¬ 
perimposed Chandra X-ray contours, 0.3-0.5 keV shown in green, and 
0.5-8 keV shown in blue. While the high energy data are located in the 
nuclear region only, the lower energy channels show extended emission 
compatible with an outflowing jet. Again, the location of the jet traced 
by X-rays corresponds to the ionized gas and radio emitting gas, and 
points towards the star forming blue clump, suggesting a possible con¬ 
nection between the star forming episode and the outflowing material. 


a FWHM sys = 286 ± 6 km s 1 , while the outflowing compo¬ 
nents have a centroid shifted by vo,om = -200 + 15 km s _I and a 
FWFLM ou , = 590 + 20 km s 1 . The projected velocity at the 10th 
percentile of the overall emission-line profile of the best fitting 
model is vio = -430 km s -1 . These high velocities confirm the 
presence of outflowing gas from the nuclear region. 

Thanks to our IFU data, we are able to spatially map 
this high velocity gas: the upper panel of Fig. [6] shows the 
contours of high velocity [OIII]/15007 emitting gas, derived by 
integrating the [OIII]/15007 emission in the spectral channels 
with v ; < -400 km s ~ 1 with respect to the bulk velocity of the 
systemic component. The contours are superimposed to the Her 
map. 

NGC 5643 is also known to host a double sided diffuse ra¬ 
dio jet, aligned with the bar, which extends by ~ 15” =1.2 kpc 
in both the E and the W direction. The radio jet is observed by 
Morris et al. (119851 1 at 6 and 20 cm, and by Leipski et al. (12006b 
at 3.5 cm. The higher resolution (1.7” x 1.7” beam) 3.5 cm VLA 
contours are shown in Fig. [6] central panel, superimposed to the 
Ha emission line map from our MUSE data. The radio emission 
morphology is similar to the EELR, although the Her emission is 
suppressed in the W direction due to the higher extinction. Given 
the presence of a low luminosity radio jet, the high velocity ion- 
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ized gas could be accelerated by a sort of “snowplough” effect, 
with the high velocity shocks induced by the outflowing plasma 
heating and compressing the external gas (see e.g. Axon et al. 
fl998l Capetti et al. H999l . 

An alternative mechanism to explain the detected radio 
emission could be the interaction between an accretion disk 
wind and the surrounding interstellar medium (see e.g. Faucher- 
Giguere & Quataert 12012, , Zakamska & Greene 120141 Nims et 
al. 120151 Harrison et al. 12015b . We note in fact that the radio 
luminosity of NGC 5643 is vL v [8.4 GHz] = 0.55 • 10 21 erg s" 1 
(Leipski et al. 120061 ), requiring a kinetic energy of the outflow 
of L wm d ~ 2 • 10 39 erg s' 1 if the efficiency of conversion in AGN 
driven winds is comparable to that in supernova driven winds. 
Thus the kinetic luminosity of the wind would be ~ 1% of 
the bolometric luminosity oif NGC 5643, making this scenario 
viable. In any case, the radio emission is consistent with being 
produced by a low luminosity radio jet or quasar wind that are 
interacting with the ISM in the host galaxy. 

We also show the X-ray images obtained filtering the Chan¬ 
dra data in the 0.3-0.5 keV and 0.5-8 keV energy bands in Fig. [6} 
bottom panel, where the high energy data are shown in blue 
and the lower energy data are shown in green. The high energy 
X-ray emission is concentrated in the central region, mostly in 
the unresolved central point source, with a minor fraction spread 
over a few arcsec around the centre. This is analogous to what 
has been observed in other similar sources (e.g. NGC 1365, 
Wang et al. [20091 NGC 4151, Wang et al. l20T0l l. The unresolved 
emission is due to the direct emission of the central AGN 
(including the primary emission and the reflected component by 
the circumnuclear gas on parsec scale). The spatially resolved 
hard emission is expected to be due to a combination of reflec¬ 
tion of the primary AGN emission by farther circumnuclear 
gas, at a distance of a few hundred pc, and the contribution of 
galactic X-ray binaries. Disentangling these components would 
require a higher S/N observation than available here, and will 
become possible with a forthcoming long Chandra observation 
scheduled in 2015. The soft emission is instead dominated by 
the extended, spatially resolved component, and is due to the 
emission of hot gas (typical temperatures l<T’~0.2-1 keV) either 
shock-heated, or photoionized by the central AGN. In particular, 
the spatial coincidence of the soft X-ray emission and the [OIH] 
emission points to a common origin of the two components, 
both due to photoionization by the central AGN of a two-phase 
gas, with the denser one associated to the [OIII] emission and 
the less dense, more photoionized one responsible for the X-ray 
emission (e.g. Bianchi et al. 120061 Balmaverde et al j20l2 ). 

The outflow direction as detected by radio. X-ray and ion¬ 
ized gas dynamics is pointing towards the location of the blue, 
Hu bright clumps, suggesting a mutual connection between the 
star formation in the clumps and the outflowing gas. In fact, we 
have seen in Sect. |2.2| that the star forming clumps A and B are 
not only on the path of the material outflowing from the nucleus, 
but they are also located at the receding edge of the dust lane 
in the bar. Bars are in fact preceded by dust lanes in their mo¬ 
tion, which have been recognized as related to shocks in the gas 
flow, and correspond to the regions of highest gas density (e.g. 
Athanassoula |1992j . Despite the high density, the gas flow is ac¬ 
companied by very strong shear that inhibit star formation. Such 
shocks are usually found in the leading edge of the bar, roughly 
parallel to its major axis. H II regions and star formation are in¬ 
stead preferentially offset from the dust lane, towards the leading 
side, as well as at the end of the bar major axis, where secondary 
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Fig. 7. Schematic view of the main structures revealed by MUSE in the 
central region of NGC 5643, drawn over the three colors image from 
Fig. □ panel c. The ionization cone borders are highlighted with red 
lines, the dust lane and nuclear dust structure location in brown. The 
two star forming clumps A and B are marked with blue circles, and the 
dashed blue ellipse, corresponding to a circle with the same inclination 
of the galaxy, shows the inner radius of the star forming ring around the 
nuclear region of the galaxy. 

gas density enhancements are seen both in simulations and ob- 
servationally (e.g. Sheth et al. 120021 . As clumps A and B are 
located on the trailing edge of the dust lane, they are not at the 
expected position for H II regions surrounding a bar, reinforcing 
the hypothesis of a mutual connection with the nuclear outflow. 
In this scenario, the outflow has been compressing the gas in 
the dust lane at the intersection point, triggering star formation 
(“positive feedback”). 

5. Conclusions 

We have presented MUSE integral field data of NGC 5643, a 
local barred Seyfert 2 galaxy at a distance of ~ 17.3 Mpc. The 
large wavelength coverage allows us to study the main optical 
emission lines in the spectra, and map the ionization and dy¬ 
namics of the gas. 

We detected a double sided ionization cone due to AGN ra¬ 
diation, almost parallel to the galaxy bar. The cone may be col¬ 
limated by a dusty structure of 580 x 125 pc around the nucleus, 
connected with the dust lane in the bar and dominated by shock 
ionization. At the center of the cone an outflow is revealed by 
the high velocity, [OIII] emitting gas, as well as in radio and X- 
ray data. The ionized gas is moving away from the nucleus at 
projected velocities up to vio = -432 km/s. 

Two Her bright, star formation dominated clumps, located at 
the receding edge of the dust lane in the bar, where the gas den¬ 
sity is highest, are located at the location where the fast outflow¬ 
ing gas from the central AGN encounter the dense material on 
the bar, strongly suggesting a relation between the two phenom¬ 
ena. The clumps are the youngest in the MUSE field of view, as 
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C„ & Evans, N. J. 2012, ARA&A, 50, 531 
Groves, B., Kauffmann, G., & Heckman, T. 2006, MNRAS, 


estimated from their Hcr(EW), and are the only star formation [2003] Kaulfmann, G., Heckman, T. M., Tremonti, C., et al. 2003, MNRAS, 346. 
dominated regions closer to the nucleus (~ 1.2 kpc) than the star 1055 
forming ring located at the edge of the field of view (~ 2.3 Kenmcutt, R. 

see Fig. |7|. We propose that the star formation in the clumps 372 961 
(SFR=0.03 M Q /yr for clump A) is due to ‘positive feedback’ of [2004] Klamer, I. J., Ekers, R. D„ Sadler, E. M„ & Hunstead, R. W. 2004, ApJ. 
the AGN outflow, that is compressing the gas at the dust lane 612, L97 
edge inducing star formation (see Fig. [7]). 

The presented data suggest that positive feedback may be 
present also at low AGN luminosity and outflow energetic. The [ 2012 ] 
upcoming observations in the framework of the MAGNUM sur-[20l4] 
vey will allow to assess if it represents an important ingredient in [2010] 
the complex black hole-host galaxy coevolution also in the local 
Universe. 
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